Capillary rogue waves 



M. ShatsQ H. Punzmann, and H. Xia 
Research School of Physics and Engineering, The Australian National University, Canberra ACT 0200, Australia 

(Dated: March 2, 2010) 

We report the first observation of extreme wave events (rogue waves) in parametrically driven 
capillary waves. Rogue waves are observed above a certain threshold in forcing. Above this threshold, 
frequency spectra broaden and develop exponential tails. For the first time we present evidence 
of strong four- wave coupling in non- linear waves (high tricoherence) , which points to modulation 
instability as the main mechanism in rogue waves. The generation of rogue waves is identified as 
the onset of a distinct tail in the probability density function of the wave heights. Their probability 
is higher than expected from the measured wave background. 



PACS numbers: 47.27.Rc, 47.55.Hd, 42.68.Bz 

Extreme wave events, also referred to as rogue or freak 
waves, are mostly known as an oceanic phenomenon re- 
sponsible for a large number of maritime disasters. These 
waves, which have heights and steepness much greater 
than expected from the sea state [1 , have become a fo- 
cus of intense research in the last decade. Recently, large 
wave events were observed in other wave systems: optical 
rogue waves [2 and acoustic waves in superfluid helium 
[3]. These discoveries indicate that rogue waves may be 
rather universal. 

The rarity of these events in the ocean and obvious dif- 
ficulties in their systematic characterization restrict the 
development of theoretical models of the rogue wave gen- 
eration. Currently four to five competing hypotheses of 
rogue waves are considered |4]. Among them is the non- 
linear mechanism of the wave- wave interactions, such as 
modulation instability. However, no direct evidence in 
support of this hypothesis has yet been found. 

Further progress in the rogue wave studies can be made 
by reinforcing that they constitute a new class of wave 
phenomena observed in different physics contexts. An ul- 
timate goal of the research into the physics of rogue waves 
is to establish mechanisms responsible for their genera- 
tion in order to (a) reliably predict their probability (in 
engineering applications), (b) generate and control such 
waves (e.g. in optical fibers), and (c) to avoid or suppress 
them (in the ocean). 

Among the key questions related to the studies of rogue 
waves are: (1) does the probability of rogue waves depend 
on the surrounding wave amplitudes, and (2) are there 
measurable characteristics of the wave background which 
would be indicative of the rogue wave probability? Mod- 
ulation instability (see [5 and references therein) has long 
been considered a likely mechanism for the generation of 
rogue waves owing to its universality. It has been found 
in surface gravity waves (Benjamin- Feir instability [6 ), 
Rossby and drift waves ([71 |8] and references therein). 
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Langmuir waves in plasma, and optical waves [5^ . Modu- 
lation instability is considered an important factor affect- 
ing the probability density function of the surface eleva- 
tion in surface gravity waves [9 . It is also a major factor 
in the generation of optical rogue waves due, for example, 
to collisions of breathers found in numerical simulations 
of nonlinear Schrodinger equation [10 . Though on the 
conceptual level there is an agreement on the importance 
of modulation instability for the rogue wave generation, 
there is no direct experimental evidence of the relation 
between the instability and probability of rogue waves. 

In this Letter we report the first observation of extreme 
wave events in nonlinear capillary waves. We show that 
the probability of capillary rogue waves strongly increases 
in a nonlinear stage of modulation instability. We present 
the first computation of the tricoherence which suggests 
that the degree of coherent four-wave coupling at this 
stage is high. 

Capillary waves belong to the higher-frequency branch 
of the surface waves, for which the restoring force is 
the surface tension. Their wavelengths are shorter than 
about 10 mm. Due to their smaller scale, capillary waves 
can be studied under well-controlled conditions in the 
laboratory using a variety of experimental methods to 
describe their space-time statistics and nonlinear wave- 
wave interactions. Capillary waves in our experiments 
are excited parametrically in a vertically shaken con- 
tainer filled with water whose depth is much greater 
than the wavelength [IT]. In these experiments waves 
are forced by shaking the container at the frequency of 
60 Hz with the acceleration in the range of a = (0.3 - 
5)^. It should be noted that the results reported here 
are not sensitive to the choice of forcing frequency. The 
strength of forcing is characterised by the value of super- 
criticality above the threshold of parametric excitation 
ath^ £ = ci/ath — 1- The strongest parametrically excited 
wave is the first subharmonic of the forcing frequency at 
/o = 30 Hz and several of its harmonics of smaller am- 
plitudes. We developed several techniques to capture the 
dynamics of nonlinear capillary waves in a broad range 
of amplitudes with high temporal and spatial resolution. 
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FIG. 1: Experimental setup: (a) Fast video recording of the 
fluorescing water surface, (b) local measurement of the laser 
beam transmission through diffusing liquid. 

In the most direct method the perturbed surface is vi- 
sualized by adding a small amount of a fluorescent dye 
(Rhodamine B) to the water and by illuminating it us- 
ing a thin (0.2 mm x 16 mm) green laser sheet covering 
almost two wavelengths at the main frequency. The mo- 
tion of the (orange) fluorescent surface in the vertical 
cross-section is then captured by fast video camera at 
600 frames per second, Fig. [TJa). An orange optical fil- 
ter in front of the camera lens discriminates against the 
direct reflection of green laser light into the camera. The 
surface contour is detected using standard edge detec- 
tion techniques to produce time-resolved records of the 
surface elevation. This technique works reliably at high 
wave amplitudes for waves of arbitrary steepness. 

Figure [2] shows a time trace of the surface elevation 
r]{t) measured at the strongest forcing, e = 5. This trace 
illustrates an extreme wave event (> 6 mm wave crest 
height). The peak amplitude exceeds the standard devi- 
ation of the wave background by a factor of more than 
five, as indicated by a horizontal grey line in Fig. |2|a). 
Two movie frames show the waveforms: before the peak. 
Fig. |2|b), and during the large event. Fig. [2|c). The 
rogue wave is characterized by an almost vertical wave 
front. Fig. [2|d) shows the probability density function 
(PDF) of the normalized wave crest heights x = ijc/cr {a 
is the standard deviation) recorded for 300 s, or 10^ wave 
periods. Up to the crest heights of x = 5, the PDF is ap- 
proximately exponential, ^ e~^^ . However the strongest 
waves, X > 5, have a probability which is substantially 
higher than expected from the e~^^ trend. 

An important question related to the high probabil- 
ity of large events seen in Fig. |2|d), is whether there 
is a threshold for the occurrence of rogue waves. The 
fast video technique does not provide sufficient spatial 
resolution to resolve wave heights at low forcing. To 
characterize the onset of rogue waves, a more sensitive 
technique is used, based on the measurement of the in- 
tensity of light transmitted through a layer of diffusing 
liquid, whereby the intensity is proportional to the sur- 
face height |T2l [13]. A thin (<0.2 mm diameter) laser 
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FIG. 2: Phenomenology of capillary rogue waves at high level 
of forcing, £ = 5. (a) Time trace of the surface elevation show- 
ing an extreme wave event in parametrically driven capillary 
waves. Video frames show waveforms (b) four wave periods 
before the large event, and (c) during the large wave event 
(gradient of the light intensity inside steep waves is due to 
the light refraction), (d) Probability density function (PDF) 
of the wave crests versus normalized crest height. 

beam is launched into the container from below, normal 
to the free surface of the water, as shown in Fig. [ijb). 
The transmitted laser light is collected onto the diffusive 
screen, which is imaged into a photomultiplier tube. This 
technique is very sensitive to small surface perturbations 
and it complements fast video imaging required at large 
wave amplitudes. 

Figure [3] shows spectra of parametrically excited waves 
at 30 Hz along with the corresponding PDFs for three lev- 
els of forcing starting at £ = 0.2, Fig.jsj^a). As the forcing 
is gradually increased we first observe the development 
of discrete frequency sidebands. Fig. |3jb), in agreement 
with classical phenomenology of modulation instability 
(see e.g. [5 ). At £ = 1.43, Fig. |3jc), the spectrum be- 
comes continuously broadened showing exponential tails. 
Further increase in forcing leads to further broadening. 
The shape of these spectra ("triangular" when plotted in 
log-linear scale) can be approximated by the hyperbolic 
secant function ^11^ . 

A new observation here is that the broadened spectra 
showing exponential tails occur above a certain forcing 
threshold (e > 1.2). Even more remarkable is the obser- 
vation of the strongly increased probability of the large 
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FIG. 3: Spectra and probability distribution functions above 
the threshold of the parametric wave excitation, (a-c) Fre- 
quency spectra around 30 Hz and corresponding (d-f) prob- 
ability density functions at four different levels of supercriti- 
cality £. 



wave events, as seen in Fig.|3[f). The onset of the tails in 
the PDF of the waves with the crest heights in excess of 
6(7 coincides with the formation of exponential spectra. 
This high probability of rogue waves is then sustained 
up to the highest levels of forcing, up to £ = 5, as in 
Fig. id). 

In the time domain the onset of discrete frequency side- 
bands corresponds to a relatively smooth wave modula- 
tion. Fig. [4]^a). Above the threshold of e > 1.2 much 
deeper modulation seen as the ensemble of individual 
wave envelopes is observed. Fig. Qb). This transition 
occurs above the same threshold as the transition to ex- 
ponential spectra. Indeed, the inverse Fourier transform 
of the hyperbolic secant spectra Ef oc sech[b{f — /o)] is 
given by s{t) = {n /b)sech[7r'^ / {bt)]e'^^^^ , which is a well 
known solution of the nonlinear Schrodinger equation, 
describing envelope solitons (see e.g. [Ml |T5] and ref- 
erences therein). Fig. shows a zoomed-in envelope 
soliton. Similar envelope solitons are also found in sur- 
face gravity waves [TBI HZ] • 

The origin of the exponential spectra has been dis- 
cussed from the point of view of envelope soliton statistics 
pT] . Stronger forcing leads to the generation of shorter 
envelope solitons and broader frequency spectra. The 
generation of nonlinearly broadened spectra in surface 
wave experiments resembles spectral broadening of the 
light pulses in photonic crystal fibers studied in the con- 
text of the supercontinuum generation [18| at its initial 
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FIG. 4: Temporal evolution of the wave height at (a) s = 
0.6, and (b) e — 1.43. The inset (c) shows zoom into one of 
the envelope solitons. 

stage. The formation of the exponential spectra in opti- 
cal fibers is attributed to either the generation of solitons 
[TOl [20] , or more recently, to the nonlinear dynamics of 
the Akhmediev breathers [21]. Our results, in particu- 
lar the transition from discrete sideband spectra to ex- 
ponential spectra. Fig. [3]^b,c), indicate that the change 
in the spectra corresponds to a transition from smooth 
quasi-periodic modulation of waves to the formation of 
ensembles of envelope solitons, as shown in Fig.|4]^c). The 
rogue wave generation probably results from a process, 
similar to the collision of breathers [T^, or it is due to 
the envelope soliton interactions observed in nonlinear 
numerical models [22] . 

In any case, modulation instability appears as a char- 
acteristic feature for such a broadening. Since four-wave 
interactions of the carrier wave with its sidebands lie 
in the heart of modulation instability, one needs to de- 
tect these interactions to positively identify the instabil- 
ity in experiments. Modulation instability is the pro- 
cess in which sidebands interact with the strong carrier 
wave simultaneously satisfying the matching rules for the 
wave numbers, /ci + /c2 = /^s + /^4, and for frequencies, 
cji + = Go's + 60^4. For example, in the spectrum of 
Fig. [sjb) one can see three peaks corresponding to the 
wave, /o = 30 Hz, and its two sidebands, /i^2 = 30 ± 0.9 
Hz. The degree of the four- wave coupling can be char- 
acterized by tricoherence [23] , or the normalized trispec- 
trum, defined as 



t^(a;i,a;2,^3) 



F1F2F3I 



+2-31 



(1) 



where Fi is the Fourier component of the surface eleva- 
tion iri{t) at the frequency uJi and Fi*^2-3 is its complex 
conjugate at the frequency cji + 6^2 — If tricoherence 
is zero, it is indicative of no coherent phase coupling be- 
tween the wave quartets, while = 1 corresponds to 
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with breaking of the modulated wave into ensembles of 
envelope solitons. These findings may also help under- 
standing the rogue wave probability and mechanisms of 
their generation in other wave systems. 

The authors are grateful to N. Akhmediev for useful 
discussions. This work was supported by the Australian 
Research Council's Discovery Projects funding scheme 
(DP0881571). 



FIG. 5: Maximum tricoherence as a function of supercriti- 
cality £. 

coherent phase coupling. Fig. |5] shows the maximum 
value of tricoherence as a function of forcing computed 
for the transmitted light intensity signals whose spectra 
are shown in Fig. [Sj The level of the tricoherence is high, 

> 0.5, at the forcing above e = 1.0 indicating strong 
phase coupling in four- wave interactions. The level of tri- 
coherence computed for the test noise signal is less than 

< 0.05. Lower tricoherence at lower forcing, e < 1.0, 
may be due to the insufficient frequency resolution in the 
tricoherence computation (A/ =1.6 Hz), such that the 
sidebands at /i^2 = 30 ± 0.9 Hz were not resolved. The 
high level of tricoherence confirms the significance of the 
underlying four- wave interactions. 

It should be noted that the dispersion relation of capil- 
lary waves allows simultaneous satisfaction of the three- 
wave matching rules. This is why three- wave interactions 
are usually considered the most likely mechanism of the 
wave- wave interactions among capillary waves. The new 
and perhaps unexpected finding in our work is that the 
dominant nonlinearity here is due to modulation insta- 
bility. High level of tricoherence indicates that 4-wave 
resonances are at work here and modulation instability 
is indeed present. Further evidence in support of modu- 
lation instability will be reported elsewhere. 

Summarizing, we have shown that extreme wave events 
occur in parametrically excited capillary waves. These 
rogue waves are very steep. Their heights exceed the 
standard deviation of the surface elevation by a factor of 
five to twelve. Rogue waves are represented in the proba- 
bility density functions as a distinct tail. The rogue wave 
probability is one to two orders of magnitude higher than 
expected from the PDF of the wave background. High 
level of phase coupling in four- wave interactions supports 
the hypothesis that modulation instability is the key in- 
gredient in the rogue wave generation. This instability 
leads to a qualitative change in the frequency spectrum 
of the waves. The onset of exponential wave spectra is 
correlated with the observation of tails in the PDF. In the 
time domain the observation of rogue waves is correlated 
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